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We solve the LO DGLAP QCD evolution equation for truncated Mellin 
moments of the nucleon nonsinglet structure function. The results are 
compared with those, obtained in the Chebyshev-polynomial approach for 
x-space solutions. Computations are performed for a wide range of the 
truncation point 10~ 5 < x < 0.9 and 1 < Q 2 < lOOGeV 2 . The agreement 
is perfect for higher moments (n > 2) and not too large xq (xq < 0.1), 
even for a small number of terms in the truncated series (M = 4). The 
accuracy of the truncated moments method increases for larger M and 
decreases very slowly with increasing Q 2 . For M = 30 the relative error in 
a case of the first moment at xq < 0.1 and Q 2 — 10GeV 2 doesn't exceed 
5% independently on the shape of the input parametrisation. This is a 
quite satisfactory result. Using the truncated moments approach one can 
avoid uncertainties from the unmeasurable x — > region and also study 
scaling violations without making any assumption on the shape of input 
parametrisation of parton distributions. Therefore the method of truncated 
moments seems to be a useful tool in further QCD analyses. 

PACS numbers: 12.38 Bx 

1. Introduction 

The DGLAP evolution is the most familiar resummation technique, 
which describes scaling violations of parton densities. Measurements of 
deep-inelastic scattering structure functions of the nucleon allow the de- 
termination of free parameters of the input parton distributions and the 
verification of so called sum rules. There exist different sum rules for unpo- 
larised and polarised structure functions which refer to the moments of the 
structure functions. From a phenomenological point of view however, QCD 
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tests based on moments Jq dxx n ~ 1 F(x,Q 2 ) is unreliable. The limit x — > 0, 
which implies that the invariant energy W 2 of the inelastic lepton-hadron 
scattering becomes infinite (W 2 = Q 2 (l/x — 1)) will never be attained ex- 
perimentally. In the theoretical approach to structure functions there are 
two ways to avoid the problem of dealing with the unphysical region x — > 0. 
The first one is to work in x-space and obtain directly the evolution of par- 
ton distributions (not of their moments). Then one has integro-differential 
equations (e.g. DGLAP one) in x and Q 2 but the integration over x goes 
for x > xq. In this case an extrapolation to the unmeasurable x — > region 
is unneeded. The second way is using evolution equations for truncated mo- 
ments of structure functions f® 2 dxx 1l ~ 1 F{x, Q 2 ) instead offer full moments. 
In the usually used method of solving QCD evolution equations, one takes 
the Mellin (full) transform of these equations and obtains analytical solu- 
tions. Then after the inverse Mellin transform (performed numerically) one 
has suitable solutions of the original equations in x-space. In this way e.g. in 
a case of DGLAP approximation, the differentio-integral equations for par- 
ton distributions q(x,Q 2 ) change after the Mellin transform into simple dif- 
ferential and diagonalised ones in the moment space n. The only problem is 
knowledge of the input parametrisation for the whole region < x < 1 what 
is necessary in the determination of the initial moments of the distribution 
functions. Using truncated moments approach one can avoid uncertainties 
from the unmeasurable x — > region and also obtain important theoretical 
results incorporating perturbative QCD effects at small x, which could be 
verified experimentally. Truncated moments of parton distributions in solv- 
ing DGLAP equations have been presented in 0. Authors have shown that 
the evolution equations for truncated moments though not diagonal can be 
solved with a quite good precision for n > 2. This is because each n-th trun- 
cated moment couples only with (n + j)-th (j > 0) truncated moments. In 
[Hj the truncated moments method has been adopted to double logarithmic 
ln 2 x resummation. There is a number of papers in which the most known 
methods for solving the Q 2 evolution equations for parton distributions have 
been reviewed (see e.g. |H],[7]). Authors compare the DGLAP framework 
for the full Mellin moments method with brute-force or Laguerre-polynomial 
approaches, used for x-space version of the evolution equation. In this paper 
we compare the solutions of LO DGLAP Q 2 evolution equations written for 
the truncated Mellin moments of the structure functions with those, ob- 
tained by using the Chebyshev-polynomial method in the x-space. In both 
these approaches we compute the truncated moments dxx n ~ l F{x,Q 2 ). 
As a test structure function F(x,Q 2 ) we take two different spin-like nons- 
inglet parton distributions. We perform the computations for a wide range 
of the truncation point 1CT 5 < x < 0.9 and 1 < Q 2 < lOOGeV 2 . In the 
next section we briefly recall an idea of the evolution equation for truncated 
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moments of parton distributions. The main topic of our paper i.e. the com- 
parison of the Chebyshev-polynomial and truncated moments techniques 
in solving the LO DGLAP evolution equation for the nonsinglet structure 
function is presented in Section 3. Finally, Section 4 contains conclusions. 



2. Truncated Mellin moments of the nonsinglet structure 
function q NS (x,t) within LO DGLAP approach. 

For (full) Mellin moments of parton distributions f(x,Q 2 ) 

i 

f(n,Q 2 ) = J dxx^f^Q 2 ) (2.1) 
o 

the DGLAP evolution equation can be solved analytically. This is because 
one obtains in the moment space n simple diagonalised differential equa- 
tions. The only problem is the knowledge of the input parametrisation for 
the whole region < x < 1, what is necessary in the determination of the 
initial moments f(n,Q 2 = Qq): 



f(n,Q 2 ) = J dxx^f^Ql). (2.2) 



Using the truncated moments approach one can avoid the uncertainties 
from the region x —* 0, which will never be attained experimentally. The 
derivation of the DGLAP equations for truncated moments of parton dis- 
tributions has been presented in The evolution equations for truncated 
moments f(xQ,n,Q 2 ) are not diagonal and therefore solving this problem 
is not so easy like in a case of the full-moments technique. Nevertheless 
this method has an advantage over other approaches, based not only on the 
cut-off for unphysical region x — * 0. The technique of truncated moments 
within DGLAP approximation enables namely to study scaling violations 
without making any assumption on the shape of the input parametrisation 
of parton distributions. While the solution of the evolution equations in 
the x-space requires knowledge of inputs f(x, Qq) with many parameters 
(fitted in detailed comparison with the data), the initial values of truncated 
moments can be obtained directly by data. Following the authors of [3], 
we have found the LO DGLAP evolution equation for the truncated at xq 
Mellin moment of the nonsinglet structure function q NS (x, Q 2 ) in a form: 



dq NS (x ,n,t) a s (t) 



2, /*»"-V"<* (2-3) 
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q N (xo,n, t) is the truncated at xq moment of the nonsinglet structure func- 
tion: 

l 

q NS (x ,n,t) = J dxx n - 1 q NS (x,t), (2.4) 

where 

* = ln-* (2.5) 
ly QCD 

and 

l 

G n (^j = J dzz^P^z). (2.6) 

xo/y 

For xq = the kernel G n (xq /y) is simply equal to the anomalous dimension 

7w(«): 

l 

7w (n) = | z^P^dz. (2.7) 
o 

Expanding the G n in Taylor series around y = 1, one has 
'x \ , , 4 ~ ~ (i + jfe_i)! 
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Truncating the above expansion at order M and using the following relation 

M k M M 

EE -EE (2-9) 

fc=0p=0 p=0fc=p 

one can find that the evolution equation (|2.3j) becomes 
dq NS (x ,n,t) a s (t) 



dt 2?r 

p=o 



and 

M 



Y J C^\x )q NS (x Q ,n + V ,t) (2.10) 



G «(7)=E^n /) (-o)?/ p , (2.H) 
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where 

+ (-+* z yi w + 2+* ifft i )1 . (2 . 12) 

Note that the evolution equations for truncated moments Q2.1U)) . (|2. 12|) are 
not diagonal but each n-th moment couples only with (n + p)-th (p > 0) 
moments. As it was shown in (3] the series of couplings to higher moments 
is convergent and furthermore the value of (n + p)-th moments decreases 
rapidly in comparison to the n-th moment. Hence one can retain from (|2.10|) 
the closed system of M + 1 equations: 

dq NS (x ,N ,t) _ a s (t) (M) NS 

jI - [ L 0,No\ x O)1 [ x O,Mo,t) 



+C^(x )g 7V5 (x ,iVo + M) + - + C^ o (xo)q" b (x ,N + M,t)] 



dq NS (x ,N + l,t) a s (t) , n {M-l) , . \s, v . . 
J t = -^-i C o,N +i( x o)l ( X 0,N + l,t) 



+^1,^+1(^0)^*^0,^0 + 2,t) + ... + Cin^o+i^o)^^^,^ + M,t)] 

^ = -^^.jVo+Afteo)? a (a? , A + M,t). (2.13) 

iVo denotes the lowest moment in calculations. The above system can be 
solved numerically like a standard coupled differential equations using the 
Runge-Kutta method. We have also found an analytical solution of (|2.13j) 
in the form: 

/ N +M 
q NS (x ,i,t) = iq NS (x ,i,t )- ^ A ik (x )q NS (x , k, t ) 
\ k=i+l 

for a, =const and 



exp(^f j (x )(t-to))+ E Ak(x )q NS (x ,k,t) (2.14) 

k=i+l 



N +M 

q NS (x ,i,t)= \q NS (x ,i,t )- A ik (x )q NS (x ,k,t ) 

k=i+l 

t\ No+M 

xexp( C/J Di M) (x )ln- + £ A ik (x )q NS (x , k, t) (2.15) 



to 



k=i+l 



G 



for the running a s . Matrix elements D\- (xq) and Aij(xo) are given in Ap- 
pendix B. For details about properties of triangular matrices like D see also 
[3]. We have made sure that the results (|2.14[) - (|2.15|) agree with the solu- 
tions obtained with the help of the Runge-Kutta method. In the forthcom- 
ing chapter we compare predictions for the truncated moments q (xq, n, t) 
obtained by solving ea. ()2.13j) with those, computed in the Chebyshev poly- 
nomial approach. 



3. Results for truncated moments of the nonsinglet structure 
function q NS (xo,n,t) within LO approximation of the 
DGLAP approach. 

We solve the system of evolution equations for truncated moments (|2.13j) 
and compare the results with predictions, obtained in the Chebyshev poly- 
nomial approach. The Chebyshev polynomials technique 10 was success- 
fully used by J.Kwiecihski in many QCD treatments e.g. [2], [HI- Using this 
method one obtains the system of linear differential equations instead of 
the original integro-differential ones. The Chebyshev expansion provides a 
robust method of discretising a continuous problem. This allows comput- 
ing the parton distributions for "not too singular" input parametrisation 
in the whole x € (0; 1) region. More detailed description of the Chebyshev 
polynomials method in the solving the QCD evolution equations is given in 
Appendix A. In this paper we use two spin-like input parametrisations of 
the parton distribution q NS {x,Q 2 ] ) at Qo = IGeV 2 , namely: 

q NS (x,Ql) = ax(l-x) 3 , (3.1) 
q NS (x,Q 2 ) = a 2 x- 0A (l-x) 2 - 5 , (3.2) 

where constants a\ and a 2 are determined by the appropriate sum rules. 
More singular at small- x input (|3.2|) incorporates the latest knowledge about 
the low-x behaviour of the polarised structure functions [5]. We start our 
analysis with a simple test, where the truncation point xq = 0. Then the 
results should be of course equal to the analytical ones: 



n (n^\\ cnqqi - n) 
' VS (''-^»-'/ VS (''-^)[^p] • (3-3) 



(Q 2 ) is the running coupling and Cf depends on the number of the quark 
flavours Nf-. 

c f = \ . (3.4) 

f 11 -IN, 

Table 1 shows the analytical values of full moments q NS (n, t) for two values 
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Q' 2 


n 


q Nb (n,Q z ) 


A 0/ 

^-Cheb A) 


axil - x) 3 


100 


1 


2.112 • 10" 1 


< 4- 10" 1 


2 


2.820 • 10~ z 


< 4 • 10 - * 


3 


7.492 • 10~ d 


< 2 • 10" 1 


4 


2.732 • 10~ d 


< 3 • 10" 1 


5 


1.204 • 10~ d 


< 7 • lO" 1 


10 


1 


2.112 • 10" 1 


< 2 • lO" 1 


2 


3.296 • lO - * 


< 2 • 10~ 2 


3 


9.556 • 10~ d 


< 5 • 10~ z 


4 


3.709 • 10~ d 


< 2 • 10" 1 


5 


1.716 • 10~ d 


< 3 • lO" 1 


a 2 x-° A (l - x) 2 - 5 


100 


1 


2.112 • 10" 1 


< 2 


2 


2.098 • 10~ z 


< 5 • wr* 


3 


5.245 • 10~ d 


< 2 • 10- 1 


4 


1.902 • 10~ d 


< 3 • 10- 1 


5 


8.502 • 10~ 4 


< 2 


10 


1 


2.112 • 10" 1 


< 9 • lO" 1 


2 


2.452 • 


< 2 • 10~ z 


3 


6.691 • 10~ d 


< 4 • 10~ z 


4 


2.583 • 10~ d 


< 9 • 10~ 2 


5 


1.212 • 10~ d 


< 2 • 10" 1 



Table 1. Test of the Chebyshev polynomial method: comparison with analytical 
results of n-th (full) moments q NS (n,Q 2 ) for different Q 2 and input functions 
q NS (x,Ql). 



of Q 2 : 10 GeV 2 and 100 GeV 2 together with the percentage errors for the 
Chebyshev results Acheb%- 

A C h b% = 1 V NS ( n ' t ) (analytical) - q NS (n,t) (Chebyshev) | _ 5 ^ 

q NS (n,t) (analytical) 

Note a good agreement of the Chebyshev solutions for q NS (n,Q 2 ) in com- 
parison to the exact analytical results. The percentage error defined in (|3.5|) 
doesn't exceed 1% in a case of the flat input (|3.1j) and 2% in a case of the 
more singular at small-x input (|3.2|) . The accuracy is better for lower Q 2 , 
when the DGLAP evolution is shorter. Using the results from Table 1, we 
expect the similar precision for the truncated moments as well. Thus we 
assume that the Chebyshev method predictions are reliable with carefully 
estimated errors: 1% for the parametrisation (|3.1|) and 2% for ()3.2|) . In 
Tables 2 and 3 we compare results for truncated at xo (0.01 and 0.1 respec- 
tively) moments, obtained from (|2.15|) (FMPR) with those, found within 
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the Chebyshev approach (Cheb.). We set again two scales of Q 2 : 10 GeV 2 
and 100 GeV 2 . 





Q' 2 


n 


q(Cheb.) 


q(FMPR) 


ai(l - x) 3 


100 


1 


1.892 • 10" 1 


2.006 • lO" 1 


2 


2.812 • 10~ z 


2.817- 10~ z 


3 


7.499 • lO - * 3 


7.491 • W~ a 


4 


2.740 • 10 -a 


2.732 • 10 -a 


5 


1.212 • lO - * 3 


1.204 • lQ~ a 


10 


1 


1.951 • 10" 1 


2.015 • 10" 1 


2 


3.289 • 10~ z 


3.293 • 10-' z 


3 


9.561 • IQ~ 6 


9.556 • 10~ a 


4 


3.714 • 10~ d 


3.709 • 10-" 


5 


1.721 • 10~ a 


1.716 • 10~ d 


a 2 x- a4 (l -x) 2 - 5 


100 


1 


1.658 • 10" 1 


1.817- 10" 1 


2 


2.082 • l()-' z 


2.090 • 10~ z 


3 


5.250 • 10~ d 


5.245 • 10-* 


4 


1.907 • 10^ 


1.902 • 10 _d 


5 


8.550 • 10~ 4 


8.502 • 10" 4 


10 


1 


1.732 • 10" 1 


1.826 • 10" 1 


2 


2.437 • 10~ z 


2.443 • 10~ z 


3 


6.692 • IQ- 6 


6.691 • 10~ d 


4 


2.585 • 10~ d 


2.583 • 10-* 


5 


1.214 • 10~ d 


1.212 • 10 _d 



Table 2. Truncated at x = 0.01 n-th moments q NS (x ,n,Q 2 ) within FMPR and 
Chebyshev approaches for different Q 2 and input functions q NS {x 1 Q 2 ] ). 



Notice a quite satisfactory agreement of the both presented methods even 
for a very small value of M (4) . The accuracy of the determination of higher 
moments is better despite the fact, that less terms (M — n) are included. 
The accuracy of the truncated moments method depends on the conver- 
gence of the expansion of G n (xQ/y), which is the truncated counterpart of 
the anomalous dimension 7 w (n). Because G n (xo/y) is expanded in powers 
of y around y = 1, the small- y region (y ~ xq) in the integral of the evolution 
equation ()2.3|) is badly reproduced. Therefore the convergence is better for 
higher moments, which have a smaller contribution from the low-y region. 
Lower moments are more sensitive to the lower limit of the integration xq 
in l)2.3j) . From the other side, for sufficiently small x$, factors x l in the 
coefficients Cp^(xo) (|2.12j) make the convergence of G n (xo/y) better. Hence 
the difference between q N (xo,n,Q 2 )FMPR and q N (xo,n,Q 2 )Cheb. is 
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Q 2 


n 


q(Cheb.) 


q(FMPR) 


ai(l - x) 3 


100 


1 


9.921 • 10~ 2 


1.236 • 10" 1 


2 


2.377 • 10~ z 


2.567- 10~ z 


3 


7.229 • 10~ d 


7.368 • 10-" 


4 


2.721 • 10~" 


2.724 • 10 _d 


5 


1.210 • 10~" 


1.203 • 10-" 


10 


1 


1.138 • 10" 1 


1.294 • 10" 1 


2 


2.883 • 10-' z 


3.011 • 10~ z 


3 


9.303 • 1Q~" 


9.401 • 10-" 


4 


3.695 • 10~" 


3.699 • 10-" 


5 


1.719 • 10~ a 


1.715 • 10~ 6 


a 2 x- a4 (l -x) 2 - 5 


100 


1 


7.112 • 10~ z 


9.086 • 10~ z 


2 


1.664 • 10~ z 


1.816 • 10~ z 


3 


5.003 • 10~ a 


5.116 • 10-" 


4 


1.890 • 10~ d 


1.895 • 10-" 


5 


8.536 • 10~ 4 


8.497 • 10~ 4 


10 


1 


8.237 • 10~ z 


9.519 • 10~ z 


2 


2.026 • 10-' z 


2.131 • 10~ z 


3 


6.446 • 10~" 


6.528 • 10-" 


4 


2.568 • 10~" 


2.572 • 10-" 


5 


1.213 • 10'" 


1.211 • 10-" 



Table 3. Truncated at xq = 0.1 n-th moments q N (xo,n, Q 2 ) within FMPR and 
Chebyshev approaches for different Q 2 and input functions q N (x, Qq). 



larger for xq = 0.1 than for xq = 0.01. Furthermore, as xq — » 1, the ac- 
cordance of q NS (xo,n,Q 2 )FMPR and q NS (xo, n, Q 2 )Cheb. becomes again 
better because of the vanishing structure functions in this limit. Compar- 
isons of q NS {x ,n,Q 2 )FMPR with q NS (x , n, Q 2 )Cheb. function of 
xq for first (n = 1) and second (n = 2) moments are shown in Figs. 1,2. 
In Figs. 3-6 we present the Q 2 dependence of q NS (xo,n,Q 2 )FMPR and 
q NS (xo, n, Q 2 )Cheb. at fixed xq = 0.01, 0.1 and for n = 1, n = 2 respectively. 
The plots are given for different M and both parametrisations (|3.1jl . (|3.2jl . 
The agreement of the truncated moment method with the Chebyshev ap- 
proach is perfect for n = 2 at xo < 0.01, independently on the inputs, Q 2 
and even value of M. The other results are also very satisfactory. The 
relative difference between q N (xo,n,Q 2 )FMPR and q (xq, n, Q 2 )Cheb. 
doesn't exceed 5% for n > 2 and not too large xo (xq < 0.1), already at 
M = 4. This difference for the first moment also decreases down to a few 
% for M = 30 and xq = 0.1 (for smaller xq the accuracy is much better). 
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0.25 




X 



Fig. 1. First truncated moment: q NS (xo, 1, Q 2 )Cheb. (solid), 

q NS (x ,l,Q 2 )FMPR (dashed M = 4, dotted M = 20) for different inputs: 
(pTTJl - upper lines and - lower lines. Q 2 = lOGeV 2 . 

0.04 | , , , — , , , — , , , — 



0.035 - 




Xo 



Fig. 2. Second truncated moment: q NS (xo, 2, Q 2 )Cheb. (solid), 

q NS (x Q ,2,Q 2 )FMPR (dashed M = 4, dotted M = 20) for different inputs: 
$FT\ - upper lines and lEOll - lower lines. Q 2 = lOGeV 2 . 



The error function 

T>M ( I q NS (xp,n, Q 2 )FMPR - q NS (x , n, Q 2 )Cheb. \ 

^ (X °' g) = q^ {xo , n ,Q2 )C heb. 10 ° % 

(3.6) 
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a 0.1 9 




1 10 100 

Q 2 



Fig. 3. First truncated at x^ — 0.01 moment: q NS (xo, 1, Q 2 )Cheb. (solid), 
q NS (x Ql l,Q 2 )FMPR (dashed M = 4, dashed-dotted M = 20, dotted M = 60). 
The upper lines correspond to the input parametrisation i|3.1[) . the lower ones to 
(ET21I . 




0.02 1 ' ' ' ' ■ — ' ' ' ' ' ' ' ' — ' ' 

1 10 100 

Q 2 



Fig. 4. Second truncated at xo — 0.01 moment: q NS (xo, 2, Q 2 )Cheb. (solid), 
q NS (x Q ,2,Q 2 )FMPR (covered with Cheb. for different M > 4). The upper line 
corresponds to the input parametrisation (|3.1f) . the lower one to l|3.2l) . 



grows very slowly with Q 2 (see Figs. 3, 5, 6). In Tables 4 and 5 we show the 
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0.14 - 




0.06 ' ' ' ' ' ' — . . . i 1 1 1 1 1 — . . . I 

1 10 100 

Q 2 

Fig. 5. First truncated at xo — 0.1 moment: q NS (xQ, 1, Q 2 )Cheb. (solid), 
q NS (xQ 1 l,Q 2 )FMPR (dashed M = 4, dashed-dotted M = 20, dotted M = 60). 
The upper lines correspond to the input parametrisation i|3.1[) . the lower ones to 
(TT21 . 



0.04 




0.015 ' ' ' ' ' ' — ■ ■ ' ' ' ' ' ' — ' ' 

1 10 100 

Q 2 



Fig. 6. Second truncated at xo = 0.1 moment: q NS (xo,2,Q 2 )Cheb. (solid), 
q NS (x ,2,Q 2 )FMPR (dashed M = 4, dashed-dotted M = 20, dotted M = 60). 
The upper lines correspond to the input parametrisation Ij3.1|l . the lower ones to 

(E2J. 
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error function R^f \xq,Q 2 ) for different M, Q 2 = lOGeV 2 and two values 
of xq : 0.01,0.1 respectively. Note that with increasing M the accuracy 



x = 0.01 


ai(l - x)' A 


a 2 x- UA (l-x) z - b 


M 


Ri 


R2 


Ri 


i?2 


4 


4 


<C 1 


6 


« 1 


10 


3 


<C 1 


5 


< 1 


20 


3 


< 1 


5 


<C 1 


30 


2 


< 1 


4 


<S 1 


60 


2 


« 1 


3 


< 1 



Table 4. The percentage error function R n = R^(x ,Q 2 ) defined in EHJ, for 
xq = 0.01 and different input functions q NS (x, Qq). Q 2 = lOGeV 2 , the values of n 
and M shown. 



x = 0.1 


a\(\ - x) 6 




M 


Ri 


R2 






4 


13 


5 


16 


5 


10 


10 


4 


11 


4 


20 


6 


3 


7 


3 


30 


4 


2 


5 


2 



Table 5. The percentage error function R n = R^ 1 (xo,Q 2 ) dehned in H3.6(l . for 
x = 0.1 and different input functions q NS (x, Qq). Q 2 — lOGeV 2 , the values of n 
and M shown. 



of the truncated moments method systematically though slowly increases. 
This improvement of the accuracy breaks however for larger M (M ~ 70 
at xo = 0.01 and M ~ 40 at xq = 0.1 ) because of increasing numerical 
errors. All presented above results concern the running coupling a s (Q 2 ). 
We have found also, that for the constant a s the error function R (|3.6|) 
grows approximately proportionally to the strength of a s : 

#4 - — • (3-7) 

Summarising, the LO DGLAP evolution of any truncated at xq < 0.1 mo- 
ment of the parton distribution can be reproduced with the satisfactory 
accuracy, where the relative error < 5%. 
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4. Summary and conclusions. 



Analysis of the QCD Q 2 evolution equations for truncated moments of 
parton distributions is very interesting both from the theoretical and exper- 
imental point of view. The truncated moments technique is complementary 
to the existing methods for solving the evolution equations, based on the full 
moments or a>space approaches. Apart from it refers directly to the physical 
values - moments (rather than to the parton distributions), what enables 
to use a wide range of deep-inelastic scattering data in terms of smaller 
number of parameters. In this way, no assumptions on the shape of parton 
distributions are needed. Dealing with truncated at xo Mellin moments: 
fx dxx n ~ 1 f{x, Q 2 ) one can also avoid uncertainty from the unmeasurable 
very small x — > region. 

In this paper we have compared the solutions of LO DGLAP Q 2 evo- 
lution equations written for the truncated Mellin moments of the structure 
functions with those, obtained by using the Chebyshev-polynomial tech- 
nique. In both these approaches we have calculated numerically and semi- 
analytically the truncated moments f£ dxx n ~ 1 F(x, Q 2 ). As a test structure 
function F(x, Q 2 ) we have taken two different spin-like nonsinglet parton 
distributions. The computations have been performed for a wide range of 
x (1(T 5 <x < 0.9) and Q 2 (1 < Q 2 < lOOGeV 2 ). Treating the Chebyshev 
results as exact, we have found that the truncated moments method is very 
promising, for any moment, together with the first one. The precision of the 
truncated moments approach is perfect for higher moments (n > 2) and not 
too large the truncation point xq (xq < 0.1), even for small M = 4. Larger 
values of M (e.g. M = 30) enables to obtain a quite satisfactory accuracy 
(the relative error < 5%) also for the first truncated moment. The original 
truncated moments technique [3] has been developed in what could im- 
prove the numerical efficiency. This technique can be a valuable tool e.g. in 
determination of the contribution to the moments of the gluon distribution 
from the experimentally accessible region. We think that the method of 
truncated moments can be useful in further theoretical and experimental 
QCD investigations. 



Chebyshev polynomial expansion within LO DGLAP evolution equations. 
In order to solve the integro-differential evolution equation 



Appendix A 




(A.l) 
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one has to expand functions q NS (x, t) into the series of the Chebyshev poly- 
nomials: 

„ JV-1JV-1 

q NS (x,t) ^q NS (x',t) = - E Viq^ix^mij^ix'), (A.2) 
where 

f 0.5 for i = , . „s 

Ui = \l for i>l ' (A - 3) 

z' = - 1. (A.4) 

in x m i n 

Ti(x) is the Chebyshev polynomial, defined as jlU| : 

Tj(x) = cos(i arccos(x)) (A. 5) 

and T}~ are nodes (zeros) of the T n \ 

2k + 1 .. , 

Tfc = cos 7T, k = 0, 1, 2, n — 1, (A. 6) 

2n 

** = *S* +1) - (A-7) 

x m j n in (|A.4|) and (|A.7|) is the smallest value of Bjorken x, involved in the 
analysis. In our computations x m i n = 10 -6 . Transformation QA.4|) converts 
the physical x-region: [x m i n ; 1] into the x' £ [— 1; 1] one, suitable for the 
Chebyshev approximation. Integration over z in the evolution equation 
(|A.1|) with q N (x,t) expanded according to (|A.2|) leads to the system of 
linear differential equations: 

= N JlH m NS (x ,t). (A.8) 

3=0 

This system can be solved by using the standard Runge-Kutta method with 
initial conditions given by the input parametrisation q NS (xj, to). N in the 
polynomial expansion ()A.2|) is equal to 20. 



Appendix B 

Analytical solution of the system of DGLAP evolution equations for 
truncated moments q NS (xo,n,t). 
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The closed system of M + 1 DGLAP evolution equations for truncated 
moments q (xo,n,t) (|2.13|) can be rewritten in the form: 

^ g y°'*> = ^[<\(x ),-(, ,iVo,t) 
+ D< & +1 (xo)q NS (x , N + l,t) + ... + D^ NQ+ M(xo)q NS (x , N + M, t)\ 



dq NS (x ,N + l,t) g a (t) ( A f_i) NS 
^ = ^^[^)v +i;wo+i( x o)g (so, iVo + 1, *) 



dq NS (x ,N + M,t) _ a a (t) (0 ) , ^5, 



^^]voWo+M(^o)r"(xo,iVo + M,t). (B.l) 

iVo denotes the lowest moment in calculations and the matrix elements 
D§\x Q ) are related to the C^\x ) (t2~T2l via 

DW(x ) = \ C S-i4 x o) 3>i . (B.2) 
J 10 j < i 



dt 



D is a triangular matrix and therefore (jB.lj) can be solved analytically using 
the diagonalising matrix A: 

j'-i 

Dij(x ) - E ^fej(a;o)^ife 

^•(x ) = n , fc 7 +1 n , ; • (B.3) 

Djj{x ) - D u (xo) 

In this way one obtains the recurrence solutions (|2.14|) . (|2.15j) . 
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